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Abstract: We investigate theoretical properties of beams of light with 
non-uniform polarization patterns. Specifically, we determine all possible 
configurations of cylindrically polarized modes (CPMs) of the electro- 
magnetic field, calculate their total angular momentum and highlight the 
subtleties of their structure. Furthermore, a hybrid spatio-polarization 
description for such modes is introduced and developed. In particular, 
two independent Poincare spheres have been introduced to represent 
simultaneously the polarization and spatial degree of freedom of CPMs. 
Possible mode-to-mode transformations accomplishable with the help of 
conventional polarization and spatial phase retarders are shown within this 
representation. Moreover, the importance of these CPMs in the quantum 
optics domain due to their classical features is highlighted. 
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1. Introduction 

States of light with spatially inhomogeneous polarization structure, such as the radially or az- 
imuthally polarized vector beams 0] [2] [3), have lately attracted major attention both in the 
theoretical and experimental area of research in optics J4). As devices for their generation 
have become commercially available, their application in optical science and engineering are 
widespread - ranging from classical to quantum optics experiments. In classical optics, when 
strongly focussed, these complex field distributions exhibit spatially separated longitudinal and 
transversal electric and magnetic field components depending on the initial state of polarization 
J5]|6l. This means, these vector beams can be sharply focussed in the center of the optical axis 
and therefore find applications in lithography, confocal microscopy 17), optical trapping and 
tweezing (8) as well as material processing [9|. Moreover, one can gain detailed insight into 
nanoscale physics due to the intrinsic strong polarization dependence of these beams iflOl . In 
the quantum regime, specially designed spatio-polarization modes can increase the coupling to 
single ions J8). Additionally, the related orbital angular momentum states help to investigate en- 
tanglement IfTTII . Very recently it has been shown that by quadrature squeezing an azimuthally 
polarized optical beam, one can generate quantum states exhibiting hybrid entanglement be- 
tween the spatial and the polarization degrees of freedom lfl2l[T3l . 

In this work we aim at establishing a proper theoretical framework for the classical optics 
description of cylindrically polarized modes (CPMs) of the electro-magnetic field. We cover 
the subtleties in terms of their description, highlight their symmetry properties and determine 
their total angular momentum. Finally, we show a suitable way of constructing a reasonable 
tool of displaying these modes on a defined pair of Poincare spheres. In particular, in Sec. 1, we 
illustrate the rotation principles underlying the intrinsic cylindrical polarization of the investi- 
gated beams. We introduce all possible configurations of these modes and calculate the Schmidt 
rank, as this is a convincing way to determine the potential inseparability of these modes. We 
complete the theoretical investigation by deducing the total angular momentum of these CPMs 



in Sec. 3. In Sec. 4, we introduce our Poincare sphere representation for the hybrid polariza- 
tion/spatial degrees of freedom characterizing CPMs. In Sec. 5, we discuss the manipulations 
of these states of light on these hybrid Poincare spheres with the help of conventional phase re- 
tarders. Furthermore, in Sec. 6, we extend our formalism to the quantum domain and illustrate 
the importance of these modes in quantum experiments. Finally, in the Appendix, we furnish 
the reader with the mathematical tools which are essential to fully cover the topics of the main 
text. 

2. Classical description of light beams with cylindrical polarization 

2.1. Introduction to rotation principles 

To mathematically introduce the radially and azimuthally polarized vector fields, let us consider 
a well collimated, monochromatic light beam with the wavelength X=2n/k propagating along 
the z-axis. The corresponding electric field can be written as 

E(r,f) = ^u(r) e -^+u*(ry^, (1) 

where % = k(ct — z) — 7i/2 lfT4l and r = xx+yy +zz. Furthermore, u(r) is a paraxial mode 
function of the form 

u(r)=/ 1 (r)x + / 2 (r)y. (2) 

The unit vectors x and y indicate the polarization in the x and y directions and f\ (r) and /2(r) 
are two independent solutions of the paraxial wave equation: 

d 2 fi,2 . d 2 h,2 . <3/i,2 ... 
<7jr ay oz 

In our investigation we are not interested in the beam propagation, and therefore we shall 
consider the longitudinal coordinate z as fixed to a value which will be omitted in subsequent 
formulas. Therefore, we define the transverse position vector x as: 

x = xx + yy. (4) 

Given a global Cartesian reference frame {x,y} on the xy-plane, it is always possible to de- 
fine a local (depending on the position vector x) polar reference frame {r(6),d(0)} with 
x = rcos0, y = rsin0, where r € [0,°°[ and 9 £ [Q,2n[. The two coordinate systems are re- 
lated through the following relation: 

f(0) = cos0x + sin0y, (5a) 
6(9) = -sin0x + cos0y. (5b) 

Since, by definition, in the {x,y} basis one has 



v I ' J VI 
then from Eq. (f5]l it follows that 

r(d)=R(d)% 6(9)=R(9)y, (7) 
where R(9) denotes the 2x2 rotation matrix which, in the {x, y} basis, can be written as: 

/cose -siney 

v ' V sin# cos0 / 



Since rotation operations are additive and commute in two-dimensional spaces, from Eq. (|7]i it 
can be deduced that for any <p € ffi the following holds: 



R((p)r(9)=R{9 + (p)x = r{9 + (p), (9) 

R(<p)d(9) = R(9 + (p)t = e{e + <p). (io) 

These equations shows that r(9) and 6(9) co-rotate with respect to the applied "external" 
rotation R((p). More generally, it is not difficult to prove that Eqs. © may be extended to two 
sets of pairs of co- and cojmfer-rotating unit vectors {r + (9) , 6 + (9)} and {r~(9),0 (9)}, 
respectively, defined as 

r ± (9)=±R(±9)x 1 6 ± (9) = ±R(±9)f 1 (11) 

where now 

R(( P )r ± (9)=r ± (9±(p), R(f)d ± (9) = S ± (0 ± f). (12) 

The two sets of pairs of unit vectors are illustrated in Fig. Q] From these equations it is clear 
that any either co-rotating u + (r, 9) or counter-rotating u~(r, 9) vector field may be written as: 

u ± (r,9)=r ± (9)a(r) + e ± (9)b(r), (13) 

where a(r) and b(r) are arbitrary functions independent of 9 but, although omitted, still depen- 
dent of z. By definition, such fields satisfy the following condition: 

u ± (r,9±(p)=R(( P )u ± (r,9), (14) 

where with /?(^>)u ± (r, 9) we had indicated the ordinary matrix-times-vector product, namely: 

R(<p)u±(r,6)=( C0&(P -™<P)( J-^V (15 ) 
vry v ' ' \ sin(p cos <p J \ y vl J 

Eq. < TT~4T > has some interesting consequences: First, we note that, since x = rcos 9 and y = 
r sin 9, we can rewrite x = x(r, 9) and y = y(r, 9), having assumed r, 9 as independent variables. 
Then, by definition, 

x(r, 9±(p) = rcos(0 ± (p) = r(cos cos (p =F sin sin 9), (16) 

and 

y(r, 9±(p) = rsin(0 ± (p) = r(sin0cos(p±cos0sin<p). (17) 

These two equations can be rewritten in compact vector notation as: 

x(r,9±cp)=R(±(p)x(r,e). (18) 

Now, since the argument of the functions (r, 9) = u* (x) on the left side of Eq. ([Pil l is exactly 
x(r, 9 ± (p), we can use the result above to rewrite Eq. (TBI as: 

u ± (R(±( P )x)=R((p)u ± (x), (19) 
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Fig. 1. Illustrating the two sets of pairs of co- and counter-rotating unit vectors 
{f + (0), 6 + (6)} and {r _ (6), (6)}, respectively. The vector f (9) is the mirror-image 
with respect to the vertical axis y of r + (6), while 6 (6) is the mirror-image with respect 
to the horizontal axis x of B (6). 



or, equivalently, 

u + (x) =R((p)u + (R{-(p)x) 7 (20a) 
vT(x)=R((p)\r(R((p)x). (20b) 

At this point, one should remember that any 2-dimensional vector field Y(x,y) transforms un- 
der a global counterclockwise active rotation by an angle <p into the new field W(x,y) whose 
functional expression is determined according to the following rule: 

W(x)=R(<p)Y(R(-<p)x). (21) 

Thus Eq. J20ab coincides exactly with the definition of a vector field which is invariant with 
respect to a global rotation by an angle <p. This means that co-rotating beams have a polarization 
pattern invariant with respect to global rotations. This result is intuitively obvious for radial and 
azimuthally polarized beams. Conversely, the corresponding Eq. d20bb for the counter-rotating 
beams has a less straightforward geometric interpretation. For this case, one can show that Eq. 
( I20bl ) can be recast in the form of a global rotation as: 

R(2(p)u-(x)=R((p)vT(R(-(p)x). (22) 

This equation simply shows that a global rotation by an angle <p upon the counter-rotating 
modes is equivalent to a local rotation by an angle 2q> on the same mode. 

2.2. Deduction of cylindrically polarized modes 

As we are now knowledgeable about all rotation properties which the cylindrically polarized 
optical beams have to obey to, we want to explicitly derive all possible cylindrical configura- 
tions of these mode functions which fulfill Eq. (TBI) taken with the "+" sign. However, it will be 
eventually shown that these coincide with the well-known radially and azimuthally polarized 
beams. 



After inverting Eq. (0 and inserting the result into Eq. (0, we obtain as a paraxial mode 
function 

u(r,0)=f(cos0/i+sin0/ 2 ) + §(-sin0/i+cos0/ 2 ). (23) 
By comparing Eq. (l23l with Eq. (l3[ , one obtains at once: 

/i(r) =a(r) cos O-b(r) sin0, (24a) 
/ 2 (r) =a(r) sin0 +b(r) cos 0. (24b) 



Since x = rcos and y = r sin 0, Eqs. ( 124al i and ( 124-bb can be rewritten as: 

/i(r)=*a(r)-y/J(r), (25a) 
/ 2 (r)=ya(r)+^j3(r), (25b) 

where Ot(r) and J3(r) can be determined by keeping in mind that f\ and / 2 need to satisfy the 
paraxial equation ([3]). In particular, we want to express the functions f\ (r) and / 2 (r) in terms 
of the first-order Hermite-Gauss ifTBI solutions of the paraxial wave equation which can be - 
apart from an irrelevant normalization factor \/%Jn - written at z = 0, as 

V/io(x)=x e -^ + > ,2 )/ w o/w , (26a) 
y m {x)=ye-^ + ^l</w Q , (26b) 

where wo is the beam waist. Thus, from Eqs. ( 125 at and d25bb and the fact that x °= y/io an d 
y oc y/ oli it follows that: 

/i(r) = -^(Avno-2»Vfti), (27a) 
/ 2 (r) = -^(B V r 10 +A W)1 ) ! (27b) 

where A, B are arbitrary numerical constants. In general, with y/ nm = y/ n m(x,y,z) (n,m € No) 
we denote the Hermite-Gauss solutions of the paraxial wave equation of the order N = n + m 
|[T5l . Consequently, all possible solutions which fulfill the global law of Eq. (TBl i with the "+" 
sign, obviously have the form of: 

u= -^=(Av/io-BV / oi)x+-^=(Bv / io+Av/oi)y- (28) 

By choosing either {A = l,B = 0} or {A = Q,B = 1} and inserting this in Eq. (1281 . one obtains 
after normalization: 



*R =-/=(Vio*+V , otf), (29a) 



1 

= -^=(-v/oix+ vioy)- (29b) 

The same steps that lead us to Eqs. ( |29at and ( |29bt could be repeated for the counter-rotating 
modes. However, it is more instructive to note that, from a visual inspection, one can immedi- 
ately deduce that and u| live in a four dimensional space spanned by the basis formed by 
the Cartesian product of the { y/jo, Vol } mode bases and the polarization vectors {x,y}: 



{Vio, Vbi}®{x,y} = {vfio£,Vio&Voit,ybi£}-- 



(30) 



However, the two vectors {u^,u|} only form a two dimensional basis and consequently two 
more vectors have to be added to span the whole space. By applying a unitary transformation 
to such a basis, we can obtain a new set of basis vectors, namely {u^ }, where: 



u R = —(-ym + mf), (31a) 



1 

flj = -^=Ooi£+Vio£)- (31b) 

The vectors {u^ , fit , ft^ , 67 }, whose intensity and polarization profiles are illustrated in Fig. [2] 
form a complete four-dimensional basis, orthonormal with respect to the scalar product defined 

as: 

2 oo 

(u(x),v(x)) = £ [ 4{x) Vi {x)dxdy, (32) 

where u(x) and v(x) are two arbitrary vector functions. It should be stressed that the basis 
vectors and do fulfill the rotation law of Eq. d 141 1. but with the "— " sign: 

ir(r,0 -a)=R(a)\T{r,d). (33) 

Furthermore, Eqs. ( TBI and d33l imply that these two sets of modes do not mix under rotation. 
This is interesting, as many recent experiments show |6][T2) that co- and counter-rotating modes 
do not occur simultaneously. 

2.3. Schmidt rank of radially and azimuthally polarized modes 

Another important property of CPMs can be characterized by the so-called Schmidt rank of 
the modes. The Schmidt rank lfl6l is an intriguing quantity as it provides information about 
potential inseparability of these modes. In the following, we demonstrate how to determine this 
rank for the radially and azimuthally polarized optical beams. This can be similarly calculated 
for the counter-rotating modes. 

If one has a look at the structure of these basis vectors, one can easily see that they are 
clearly non-separable, in the sense that it is not possible to write any of these vector functions 
as the product of a uniformly polarized vector field times a scalar function. Such a "degree of 
non-separability" may be quantified by the Schmidt-rank of the function, remembering that all 
separable functions have a Schmidt rank of 1 . Conversely, we are now going to prove that 
and u| admit a Schmidt decomposition with rank 2: 

Let ej and be two orthogonal unit vectors defined as ei = x and §2 = y, with 

e?-e y = %, j,ye{l,2}, 
and let vi and v 2 be two orthogonal functions defined as 

vi = Vio(x, y,z), v 2 = \l/ 0i (x,y,z), (34) 
where, because of the orthogonality of the Hermite-Gauss modes, one has: 

v*vj dxdy = € { 1 , 2}. 

Then we can immediately write down, for example for the mode: 

2 



E(A,-) I/2 e,-v„ (35) 



= 1 



(a) 



(b) 
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Fig. 2. Complex polarization patterns of (a) u^, (b) u|, (c) u s , (d) u A , underlayed with 
the doughnut shaped intensity distribution. 

with Ai = A2 = 1/2. The right side of Eq. (f35T > represents the Schmidt decomposition of the 
vector function u^ with rank 



If =2. 



For the function u| the same demonstration holds but with vi and V2 defined as 

vi = -yoi(x,y,z), V9 = \jf l0 (x,y,z), 

This completes the proof. 

More generally, the substitution of Eq. ( I29ab or ( I29bl i respectively into Eq. (fl4l) gives 



(36) 



(37) 



*ab = — iVw (Ax + Bf) + y/bi (— Bx+Ay)] 



1 

1 e 1 * 

-^=fiy/io + -^=f 2 ybi, 



(38) 



where we have chosen |A| 2 + \B\ 2 = 1 to guarantee normalization: (uab,uab) = 1. From Eq. 



it follows that u^ B has Schmidt rank 2 whenever 

fl-f 2 =AB*-A*B = Q, 



(39) 



namely when lm(AB*) = 0. Conversely, u^ B becomes separable when f 2 = Ct\, namely when 
B = CA, A = CB, where C is a proportionality constant. The substitution of the latter equation 
in the previous one furnishes — B = C 2 B which implies C = ±;, namely A = ±iB. In this case 
it is easy to see that Eq. (f36t reduces to 



(40) 



which means that the only separable cylindrically symmetric co-rotating modes are the circu- 
larly polarized Laguerre-Gauss modes with one unit (\£\ = 1) of orbital angular momentum: 
= ( V10 i 'Vol )/v2~- It should be noticed that it is only thanks to the opposite double sign 
in the last two terms of Eq. (l40l > that condition (fl4] > is satisfied. 



3. Angular momentum of cylindrically polarized optical beams 

To complete our discussion about the properties of the "+" and "— " vector bases, we note that 
all four vector functions {u^",u^,u^,u^} have the total (orbital plus spin) angular momen- 
tum (TAM) of zero which is deduced below. Following Berry |[T7l . we can write the "local" 



(position-dependent) density of the linear momentum carried by the paraxial mode u of the 
electromagnetic field as: 



p(r) = Im[u* • (Vj.)u] + ilm[Vj_ x (u* x u)] 

= Porb(r)+p sp (r), (41) 

where Vj_ = kd x +fdy and u* ■ (Vjju = £j«*(Vj_Mi). A prefactor of c 2 £o/(2co) has been omit- 
ted. The symbol "x" denotes the ordinary cross product in M 3 . In Eq. d4l"l ). the first term p rt>( r ) 
is called the orbital Poynting current and is independent of the polarization state of the beam. 
Conversely the second term, the so-called spin Poynting current p sp (r), contains the local "den- 
sity of spin" of the beam via the term Im(u* x u) ifTTl . An explicit calculation from Eq. (|4"TT i 
furnishes: 

Px (r) = ilml/iWO +fi(d x f 2 )] + - k \m[d y {flf 2 )] , (42a) 

Py (r) = ilm [f x (d y fi ) + f x (d y f 2 )] - ilm [d x (ftf 2 )\ , (42b) 

p z (r)= \fi\ 2 + \f 2 \ 2 , (42c) 

where k indicates the wave vector k = co/c. The power of a light beam as actually measured by 
a photo-detector in the laboratory equals the flow of the Poynting vector of the beam across the 
detector surface. Since the electromagnetic linear momentum density p is equal to 1/c 2 times 
the Poynting vector, it follows that a quantity experimentally observable is P • z, where 



P 



J J p(r) dxdy 

Porb + Psp, (43) 



is the linear momentum of the field per unit length. Without loss of generality, we have assumed 
that the detector surface coincides with the xy-plane normal to the beam propagation axis z. The 
double integration in Eq. (l43l is extended to the whole xy-plane. This means that for any pair 
of functions f\ and f 2 , normalizable in the xy-plane, which drop to zero at infinity, the "spin" 
surface terms d y (ftf 2 ) and d x (ftf 2 ) on the right sides of Eq. (142 at and Eq. (I42bt respectively, 
do not contribute to the measured momentum. In other words one has P sp = 0. It is worth noting 
that although f\ and f 2 are both functions of r, as opposed to x we can practically evaluate Eq. 
( |43| > at z = because the quantity P is independent of z lfT8l . 
The total angular momentum density is 

j(r) = r xp(r) 

= r x [p orb (r)+p sp (r)] 

= l(r)+s(r), (44) 

where the orbital angular momentum density is 

l*(r) =y(\fi\ 2 + \f 2 \ 2 )-^lm(ftd y f l +f 2 t d y f 2 ), (45a) 

l y (r) - | 2 + |/ 2 | 2 ) + ^Im(/f d x fi +f 2 d x f 2 ) , (45b) 

l z (r) = Im [xiftdyfi +f 2 *d y f 2 )-y(ftd x f 1 + f 2 d x f 2 )] , (45c) 



and the spin angular momentum density is 



s x (r)= -|lm[^(/r/ 2 )], (46a) 

s y (r) = ^m[dy(ftf2)], (46b) 

s z (r) = - -j-Im [xd x {ftf 2 ) +yd y {fth)\ . (46c) 

We proceed as in the linear momentum case and we write the angular momentum of the field 
per unit length as 



J 



Jjj(r)dxdy 

L + S, (47) 



were, similarly to P, also L and S are independent of z lfl8l . Again, for any pair of functions 
/i (r) and /2(r) that vanish sufficiently fast for \xy\ — > °°, integration by parts leads to S x = = 
S y , and to 

JJ[xd x {ftfi) +ydy(ftf 2 )]dxdy = -2 jjftf 2 dxdy. (48) 
From this relation it follows straightforwardly that 

{ftfi-hffidxdy. (49) 



Eq. (|49b clearly reproduces the well-know expression of the helicity S z = i(aj5* — a*f5)/k of 
a uniformly-polarized beam when one chooses f\ = ai//(r) and f 2 = y3 i//(r), where is a 
normalizable function in the xy-plane and a,/3 a pair of complex number such that |a| 2 + |j3| 2 = 
1. 

At this point, we can easily evaluate J for the set of our four cylindrically polarized modes 
without performing any actual calculations. In fact, since J is independent of z, we can evaluate 
it at z = where both y/\()(x,y,z = 0) and y/oi (x,y,z = 0) are real-valued functions. Then from 
Eqs. ( 145H46I ) it follows at once that S = and that 

ydfrf + lfil^dxdy, (50a) 

Ly =- JJ jc(|/i | 2 + \f 2 \ 2 )dxdy, (50b) 
L z = 0. (50c) 

However, it is easy to see from Eqs. ( I26al i and ( I26bb that by definition | \jfio(x,y,0)\ 2 and 
IV^oit^^O)! 2 are even functions with respect to both x and y coordinates. This automatically 
implies, because of symmetry, that L x = = L y . 

We have thus demonstrated that the TAM of our modes is strictly zero. Alternatively, the 
same result could have been found in a more intuitive but less formal way by writing the modes 
{u^,u^ ,u|,u^} in the left/right circular polarization basis and in the <j>^ Laguerre-Gauss 
spatial basis. In this case, both the polarization and spatial modes carry a unit of angular mo- 
mentum which may be summed or subtracted from each other, thus leading to a straightforward 
physical interpretation. 



Fig. 3. Poincare sphere representation for an arbitrary two-dimensional system. Here, the 
radius r of the sphere is fixed to one, and spherical angles { ■&, (j)} are related to the Stokes 
parameters via the relation: {cos#,sin#cos</>,sin#sin0} = {Si/Sq,S2/So,S^/So}, where 
So denotes the total intensity of the beam. 



4. Hybrid Poincare sphere representation 

In general, fully polarized optical beams, similarly to quantum pure states of two-level systems 
lfl9l . can be described by points on a sphere of unit radius H20I . the so-called polarization 
Poincare sphere [PPS, Fig.[4](b)]. The North and South Poles of this sphere represent left-and 
right-handed polarization states that correspond to ±1 spin angular momentum eigenstates, 
respectively. Similarly, the spatial-mode Poincare sphere [SMPS, Fig. [4] (a)] 1031 sticks to the 
same superposition principle but in terms of the spatial distribution of the optical beams. In 
this case, the North and the South Pole of the sphere represent optical beams in the Laguerre- 
Gaussian modes (j)^{ , which correspond to ±1 orbital angular momentum eigenstates lED 
respectively. Each point on the Poincare sphere, both on the polarization and the spatial-mode 
one, may be put in one-to-one correspondence with a unit three-dimensional real vector S = 
{Si/So,S2/So 7 S?,/So}, where Sq = S\ + S\ + S\, as shown in Fig. [3] The three components Si, 
i G {1,2,3} are known as the Stokes parameters |j20l and Sq denotes the total intensity of the 
beam. 




Fig. 4. Fundamental idea behind the new Poincare sphere representation: Combining the 
SMPS (a) and the PPS (b) to the hybrid Poincare spheres (HPSs) (c) and (d). 



In this section, we present a novel representation for cylindrically polarized optical beams 
based on the simultaneous use of both the PPS and the SMPS representation. We show that 
it is possible, under certain assumptions, to connect the two polarization and spatial represen- 
tations to obtain two hybrid Poincare spheres (HPSs) that embody all characteristics of these 
complex modes of the electromagnetic field. Up to now, only uniform polarization patterns 



could have been displayed as points on the PPS. However, with this new representation, non- 
uniform polarization patterns are represented for the first time. This is especially helpful - let 
it be in theoretical considerations or in experiments - if one investigates effects deriving from 
these spatially dependent, complex polarization patterns. Furthermore, manipulations of states 
within this representation in accordance with the conventional tools such as the PPS and SMPS 
can be carried out. This is investigated in Sec. 5. 

In principle, the two-dimensional polarization and spatial spaces, respectively spanned by 
the bases {x,y} and { V^io, W()\}, may be combined in several different manners to produce a 
four-dimensional mixed polarization-spatial mode space. A complete representation of such a 
space would require the simultaneous display of eight real numbers which is impossible in a 
three dimensional space. Consequently, our aim is to find an appropriate way for describing 
this space by using conventional tools, such as the PPS or the SMPS respectively. We begin by 
noticing that if we are permitted to discard some information, then a lower-dimension visual 
representation may be possible. The characteristics of this picture would then be fixed by the 
kind of degree of freedom (DOF) we choose to represent. Driven by certain recent experimental 
results lfl2l . here we look for a hybrid representation, which would be suitable to describe phys- 
ical situations where co- and counter-rotating modes do not mix with each other. This absence 
of mixing occurs, for example, when the light propagates through cylindrically symmetric opti- 
cal systems as fibres, non-astigmatic lenses, et cetera. In this case, the disregarded information 
would amount to the relative phase and amplitude between the sets of co- and counterrotating 
modes. The representation deriving from these constraints is illustrated in Fig. 2] The North and 
South Pole of sphere (c) represent radially and azimuthally polarized states whose total angular 
momentum (orbital plus spin) is equal to zero. Analogously, sphere (d) represent counter-radial 
and counter-azimuthal states, still with total zero angular momentum. Each of these two pairs 
of states represent a novel kind of hybrid DOF of the electromagnetic field, since such states 
describe neither purely polarized nor purely spatial modes of the field. 

Mathematically, the combination of the SMPS and the PPS is realized by forming the four 
dimensional basis given by the Cartesian product of the the spatial mode { y/io, Voi } and polar- 
ization {x,y} bases: 

{Wo,W)i}<8>{x ) y} = {v/i x,v/ioy, roix, roiy}- (51) 

To span this four dimensional space, we choose the basis {u^,u^,u^,u^}. As it was shown 
above, these two sets of vector functions, namely the "+" and "— " ones, do not mix under 
rotation and, moreover, obey the two fundamental rotation laws in Eqs. ([Pil l and ( f33b . Such a 
physical constraint may be mathematically implemented by a superselection rule |[T6l which 
forbids interference between states separated according to the law in Eq. ( TBI . Therefore, we 
use this law to split the four dimensional space { i/io x , V'loy, Wq\%, WoiS} into the Cartesian sum 
of two subspaces spanned by the "+" and "— " sets of modes: 

{ Vio , Voi } ® {% y } = {i4 > U A } © i U R > U A } ■ ( 52 ) 

This shows, that depending upon their behavior under rotation, cylindrically polarized optical 
beams may be subdivided in two independent sets which can be represented by points on the 
surface of two distinct hybrid Poincare spheres, as illustrated in Fig.@] It is worth to stress once 
again that a beam represented by an arbitrary superposition of either the standard basis vectors 
{ Vio x , Ywf: Vol*: Yoif} or the cylindrical basis vectors {u^,u^,u^,u^} needs four complex 
numbers to be described completely. This amounts to eight real numbers which can be reduced 
to seven because of normalization. However, our superselection rule defined above reduces the 
number of real parameters necessary to describe cylindrically polarized beams further to 3 © 3, 
thus permitting the introduction of a "two-Poincare Cartesian sphere" representation. 
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Fig. 5. Polarization states on the HPSs, represented on the sphere of the "- 
the "— " modes (b). 



1 modes (a) and 



Besides evident visualization properties, our representation permits to describe straight for- 
wardly cylindrically polarized states of light, in terms of hybrid Stokes parameters, as opposed 
to the either polarization or spatial Stokes parameters. These hybrid Stokes parameters convey 
information about both polarization and spatial-mode degrees of freedom simultaneously. They 
are naturally defined as: 
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where the symbol * denotes the complex conjugation. f£ = (irj , E) and f% = (u^ , E) are the 
field amplitudes of the electric field vector in the bases {u^,u^}, where 

E = / A ± u|+/ R ± u± (54) 

These amplitudes can be expressed in terms of the spherical coordinates on the two independent 
hybrid Poincare spheres (HPSs) as 

jf= cos (0/2), (55a) 
=exp(/0)sin(0/2). (55b) 

The first sphere of the HPSs represents the "+" modes [Fig. [5] (a)] and the second the "— " 
[Fig- 0(b)] modes. It is possible to show that these Stokes parameters can be actually measured 
by means of conventional optical elements. This characteristic is particularly relevant for the 
possible quantum applications of our formalism, where simultaneous measurability of Stokes 
parameters describing spatially separated optical beams, is crucial l22l . 

Concerning the structure of the HPSs, an interesting feature is that every point on the merid- 
ian between the North and the South Pole (S* e [0,7f], ± = 0) describes a locally linear 
polarization state as illustrated in both Figs. [5] (a) and (b). Except for these points and the two 
ones on the S3 axis (9 = f,0 ~ j7i}) which represent circularly polarized states, all re- 
maining points on the spheres describe states of non-uniform elliptical polarization. This is 
fully consistent with the structure of the conventional PPS or SMPS representation. 
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Fig. 6. Manipulation of states on the HPSs, represented on the sphere of the "+" modes (a) 
and the "— " modes (b). 



5. Manipulating states on the HPSs 

As in terms of the conventional PPS and SMPS representation, conversions of states by means 
of polarization and spatial phase retarders such as waveplates or cylindrical lenses, respectively, 
are possible. In general, waveplates only act on the polarization part of the vector beam and 
cylindrical lenses on the spatial mode part, as it was shown in lEH . In terms of the PPS, two in- 
dependently adjustable quaterwaveplates (QWPs) and one halfwaveplate (HWP) are sufficient 
to control polarization l20l . so reaching every point on this sphere is possible. In analogy to 
this, manipulating states on the SMPS is possible, as well. However, one cannot use waveplates 
but instead cylindrical lenses (151 to control the spatial mode. Displaying the information from 
both the PPS and the SMPS simultaneously with the new representation, combinations of both 
waveplates and cylindrical lenses are necessary. The process of converting a state on the HPSs 
to a different one - for instance, going from Sf via Sf, —Sf and — Sf back to Sf - can be 
achieved by applying two halfwaveplates and rotating one of these as illustrated in Fig. [6] The 
fixed HWP induces a local flip of the vectors in the polarization pattern and the other one rotates 
continuously these vectors. However, due to the special cylindrical symmetry of our states, no 
"arbitrary-to-arbitrary" l23l transformations with the help of waveplates or cylindrical lenses 
are permitted. Therefore, only a few transformations can be accomplished without altering the 
symmetry of the state. These are listed below: 
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where 9 £ [0,7l] and <p refers to the angles parametrizing a sphere, as shown in Fig. [7] These 
constraints are due to the fact that phase retarders imprint a spatially uniform phase shift onto 
the cross section of the beam, which is not, per definition, cylindrically invariant. These con- 
straints are explicitly deduced in Appendix. 

To complete the discussion about possible manipulations, we notice that the two sets 
{u^,u^} and {u^,u^} are connected to each other by a mirror image symmetry operation, 
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Fig. 7. Transformations on the "+" HPS: Without altering the rotational symmetry of the 
state, these transformations are possible with the help of conventional phase retarders. 



which is illustrated in Fig. [8] This can be practically performed by a halfwaveplate l23l . So, 
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Fig. 8. By a mirror image symmetry operation one can pass from one sphere to the respec- 
tive other sphere. 



jumps between the "+" and the "— " spheres can be carried out easily and they considerably 
enlarge the number of possible transformation. 



6. Quantum properties of cylindrically polarized states of light 



As we have discussed only classical properties of cylindrically polarized states of light so far, 
we now cover their quantum aspects. In this section, we give a few examples of application of 
our hybrid polarization/spatial formalism to some quantum states of light. 

To begin with, let us write the electric -field quantum operator for paraxial beams of light in 
the same units of Eq. (Q]i as: 

E(r,f)=E+(r,r)+E-(r,f), (57) 

where E~ = (E + ) t , and 

1 2 

E + (r,f) = -X LM;UmVW ' x , (58) 

L X=\n,m 

with i\ — x and £2 = f and n,m€ {0,1,2,...}. The paraxial annihilation operators a^ nm satisfy 
the canonical commutation relations 

[ & Xtmi&\i n > m i\ = 5a A '<>„„'(>,„,„/. (59) 

In order to have a correct description of quantum noise in light beams with cylindrical polar- 
ization, we need to establish first a correspondence between the classical and the quantum rep- 
resentation of these beams. For sake of simplicity, here we restrict our attention to co-rotating 
modes solely, namely u| and u^. Having in mind this goal, let u^lr) be the generic cylindrical 
mode (the superscript "+" will be omitted from now on) that now we write as 

u AB {r) = -j=[Z(A\i/io-By i)+f(BVio+Ayoi)] 

= Aur + B\x A , (60) 

where Eqs. ( |29at and d29b| > were used in the second line and again \A\ 2 + \B\ 2 — 1 is assumed. 
Since (u^u^) = 0, then two different modes Uab(i") and u A /gi(r) satisfy the relation 

(u AB (r),u A , B ,(r)) =A*A' + B*B'. (61) 

Now the question is: What is the correct form of the field operator a^B that annihilates a photon 
in the mode U4fi(r)? We seek an answer to this question by requiring that: 

1. The coherent state \a) = exp(ad\ B — a*«Afi) |0) must produce a coherent signal ,5? = 
(a|E|a) equal to the classical electric field: 

& = \ [au AB (r)e- ix + a*u* AB (r)e ix ] . (62) 

2. The single-photon state |1) = fl^O) must generate a photon wave-function (0|E + |1) 
equal to: 

(0|E+|l) = iu AB (r) e -^. (63) 
With the help of Eq. ( |60T l we guess the following form for the sought annihilation operator a AB : 
a AB = A {-^^) +B ( ^ - ) , (64) 



where we have used the polarization indexes x for A = 1 and y for X =2. From Eqs. ( 159164b it 
follows that: 

[&a>b'^Ib] =A*A'+B*B', (65) 

which, together with Eq. doTl i. shows that a A i B i and a\ B commute whenever the two modes 
UAs(r) andu A / B /(r) are orthogonal. 

6.1. Coherent states 

Let Z),(/3), i € {1,2, ... ,°o} be the coherent-state displacement operator defined as 

A(j3)=exp(j3fl,-j3*4) ; ( 66 ) 

where the single label i embodies X,n,m, so that = L;iL„L m , an d we choose the first 4 
values of the index i in such a way that 

«1 = fl*10i «2 = A>01j «3 = &xOl, «4 = a.vlO- (67) 

If we rewrite the field mode functions as 

v,(r,f) = e A v/ nm exp(-/^), (68) 
then Eq. (|58"I ) takes the simpler form 

E+ = i£v^. (69) 

2 i=l 

Since any coherent state |j3); = D,-Q3)|0) satisfies the relation 

&-\P)i = livjaj\P)i = UiP\p)u (70) 

z ; . =1 z 

it immediately follows that if we define |cki, 0^,0^,0^) = At(a 4 )D3(a3)^(a2)I>i(ai)|0), 
then 

E + |a 1 ,a 2 ,«3,« 4 ) =i^v i 4|ai,a 2 ,a 3 ,a 4 ) 
2 i=i 

1 4 

= -£v,ff i |a 1 ,a 2l «3,a(), (71) 

which trivially implies that 

^ = (ai,a2,a3,a 4 |E + + E"|ai,a 2 ,a3,«4) 
1 4 

-£(v,a, 1 v*a,') 
z 1=1 

= - [ft (ai 1^10 + a3 Vbi) + f (<Wio + <Woi )] <?~'* 



(72) 



Comparison of Eq. ( TTZb with Eq. (l60b shows that if we choose 
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then we have: 



= exp 
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exp[aa' AB -a a AB ) 



(74) 



and Eq. (f62t is automatically satisfied. 



6.2. Fock states 

The validity of Eq. ( f63T > can be checked by defining the single-photon ladder operator A (| ) as: 



A(§) = £^-, with £|£| 2 = 1, 

1=1 1=1 

that satisfies the canonical commutation relations 

[a(£),a^)] = i. 

Let 1 1) be the single-photon state defined as |1) = A T (^)|0). Then 



(75) 



|fi + |l)=^(0|£v,-a ( -f§ y aj|0) 
i=i /=i 

= i£^Vi(0|^4|0) 

1 4 



Now, if we choose 



then we obtain 
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(76) 



(77) 



(78) 



in agreement with Eq. d63l l. Thus, we have proved that Eq. ( |64l furnishes the correct expression 
for d AB . 

6.3. Squeezed states 

As a last example, here we write explicitly the expression for an azimuthally polarized squeezed 
state. The extension to the other cylindrically polarized states is straightforward and, therefore, 
will be omitted. As a result of such calculation, we shall find that classical inseparability auto- 
matically lead to quantum entanglement when the azimuthally polarized beam is prepared in a 
nonclassical state lfl"2l . 

The "natural" definition of squeezing operator for the azimuthally polarized field is 



(79) 



where £ = se'® is the complex squeezing parameter. By substituting Eq. (l64l evaluated for 
A = 0, B = 1 into Eq. {79]), it is not difficult to obtain 



5 A (C)=5 3 (C/2)5 4 (C/2)534(-C/2), 



(80) 



where we have defined the one- and two-mode squeezing operator 



HQ 
SijiQ 



exp ic*(o/) 2 -^C(«D 2 



(81) 



exp ( C,*&icij — C^i^) 



(82) 



Thus, the squeezed-vacuum state 



5 3 (C/2)5 4 (C/2)534(C)|0>, 



(83) 



is clearly entangled, since 



534(O|0) 



1 



£( e - w tanhs)» 3 |«)4, 



coshs 



(84) 



n=0 



represents the well-known entangled two-mode squeezed vacuum Il24l . Thus, we have shown 
that classical inseparability naturally lead to quantum entanglement. 

7. Conclusion 

In conclusion, we have established a proper theoretical framework for a classical description 
of cylindrically polarized states of the electromagnetic field. We have highlighted the rotation 
principles of these modes which can be derived from their very particular intrinsic symmetries. 
Moreover, we found that all CPMs have a zero total angular momentum and are clearly non- 
separable which can be proven by their Schmidt rank. This subtlety leads to intriguing features 
of these modes in both the classical and quantum domain of optics, the latter being recently 
demonstrated in experiments fl2|. We have exploited these findings to present a new way of 
visualizing CPMs on a pair of hybrid Poincare spheres. This permits, at a glance, to easily 
recognize how phase retarders such as waveplates or cylindrical lenses act on these complex 
polarization patterns. Furthermore, we have extended our discussion of the properties of these 
modes to the quantum domain and have shown that the classical inseparability property leads 
to the intriguing feature of entanglement. 

In general, our results presented here give insight into the theoretical structure of these modes 
and their particular behaviour. We believe that the presented formalism and results are of utility 
to both the quantum and the classical communities especially in view of the recently growing 
interest in the theory and applications of light beams with complex spatial and polarization 
patterns. 

The authors thank Peter Banzer for fruitful discussions. 

8. Appendix 

In our hybrid Poincare sphere representation, only certain transformations are permitted. Due 
to the fact that these constraints are not intuitive, we demonstrate here how they can be derived. 
Generally, besides the particular cases in illustrated Fig. [7] the fundamental rules that an optical 
transformation must obey in order to not violate the global rotation law in Eq. (TT4T > and con- 
sequently not to generate mixing between the co- and counter-rotating modes, can be deduced 
from perfectly general principles, as follows. 



Consider a generic transformation & = J% x & that acts upon both spatial (^#) and po- 
larization (3?) degrees of freedom of the field separately. Here, it is convenient to adopt a 
quantum-like notation, and write 

4=|«i), ?=|«2>, (85a) 

i/a 10 (x) = |^i(x)), VoiW = (85b) 

where the symbol "=" stands for: "is represented by". Note, that the quantum-mechanical no- 
tation is just a convenient manner of writing our classical modes of the electro-magnetic field 
and it has nothing to do with quantum-mechanical interpretation. However, for the sake of com- 
pleteness, our hybrid classical description will be extended to the quantum domain in Appendix 
II. Thus, for example, our set of cylindrically polarized modes may be represented as: 

uf = |4) = -^(±|vi>|«i) + |%)|«2», (86) 

v^ = \u^) = -^=(T\Y2)\ei) + \Wi)\e2))- (87) 

In passing, we note that these states have the formal structure of the quantum Bell-states. 
More specifically, if one establishes the following formal equivalence between the two qubits 
(A,B) Hilbert space and our two degrees of freedom (polarization and spatial mode) space: 
\0a) ~ x, \ l A ) ~ y, |0b) ~ y/io, |1b) ~ Vbi> men it is not difficult to show that our basis is 
mathematically equivalent to the quantum Bell basis J25): 

{u+ u+,u R , u A } ~ ,-!*"), |¥+)}, (88) 

where |4>±) = (|0 A ,0 B ) ± |1 A , l a » /Viand 1^) = (|0 A , U) ± |U,0 B » /y/2. This analogy has 
been recently exploited in 11261 . 

Now, let us rewrite Eqs. J20ab and d20bl i in quantum-like notation as: 

|« ± (x))=G ± (( P )|« ± (x)), (89) 

where the two operators (j(<p) are determined by their action upon the generic basis state 
|y/,(x))|ej). According to Eqs. (fT~9b one has: 

G ± (( P )\Wi(x))\e J ) = \Yi(R(T<p)x))&((p)\ej}, (90) 

where the rotation operator R((p) is expressed in terms of the matrix elements Rkj((p) of the 
rotation matrix R((p) as: 

2 

%)l«y> = L%(p)ht>, (91) 

k=l 

From Eqs. ( I26ab and ( I26bb it follows that y/io and y/oi transform as the components of the 
position vector x, namely: 

i/Aio(#(T<P)x) = yio(x)cos<p± y/ 01 (x)sin(p, (92a) 
y/bi(/?(=F<P)x)= T Vio(x)sin<p + y/ i(x)cos<p. (92b) 
These formulas may be put straightforwardly in the operator form which reads: 

\Yi(R(T(p)x)) =M(T(p)\Yi(x)) 

2 

= £tf a (TP)|V;(*))- (93) 

i=i 



Note the different order of the indices in the matrix elements appearing in Eq. d9"TT l and Eq. 
d93l >. Finally, gathering all these formulas together we obtain: 

G^((p)\Y i (x))\ej)=R(T(p)\Yi(x))R(( P )\e J ) 
= Y, R u(T<P)\Yi(x))Rkj((p)\ek), 

k.l 

= Y J R i i(T<p)R] k (<P)\Vi(x))\ek), 

k.l 

= £[R(T9)®*(-9)]jy ft lwW>l«*>, 

k.l 

= R(T(p)®&(-cp)Mx))\ej), (94) 

where A T indicates the transpose of the generic matrix A, and we used the fact that for orthog- 
onal matrices one has R T ((p) = R(—<p). Here, in the direct product J?(=F<P) ® R(—<p) the first 
operator acts upon the spatial degrees of freedom solely, while the the second one acts on the 
polarization degrees of freedom. Since in all previous calculations the basis state |i//,(x))|e 7 ) 
was arbitrarily chosen, the relation d94l i above is perfectly general and may be rewritten as: 

G±((p)=R~(T(p)®R(-(p). (95) 

Now, let M and P be such that T = M ®P, the two spatial and polarization operators cor- 
responding to the physical transformation & = ./# x Then, by definition such operation 
transforms |m ± (x)) into T\u ± (x)). It is clear that this new state maintains the same symmetry 
under rotation of the original state \u (x)) if and only if: 

G ± (f)f\u ± (x)} = f\u ± (x)}. (96) 
By multiplying both sides of Eq. d89l by T, one straightforwardly obtains: 

f\u ± (x))=fG ± (< P )\u ± (x)). (97) 
Thus, by using Eq. < f97T > into Eq. d93l l one easily obtain: 

G ± (( P )f\ u ± (x))=fG ± ( ( p)\ U ± (x)). (98) 

or, equivalently: 

[G ± (( P )f-fG ± ((p)]\ U ± (x))=0. (99) 
This equation can be easily recast in commutator form by rewriting it as follows: 

[G ± ((p),f]\ u ± (x))=0. (100) 

Of course, any operator f that commutes with G ± (^>), namely [G ± (^>), f] = 0, automatically 
fulfills Eq. (1 1 00b . However, this condition is sufficient but not necessary because Eq. dlOOb 
represents a much weaker constraint, as it only requires that the kernel (or, at least, part of it) 
of the operator [G ± ((p),7'] coincides with the subspace spanned by either |m + (x)) or |m~(x)). 
In other words, all what is required to t by Eq. ( 1 100b , is that the commutator [(j(p), t ] has a 
null eigenvalue in correspondence of the eigenvector |m ± (x)). 

Now, having this caveat in mind, we can rewrite the condition ( 1 1001 ) as follows: 



= [M®P] [R(T(p)®R(-(p)], 



(101) 



or, equivalently: 



R( T( p)M[R( T (p)] ®R(-(p)P[R(-(p)] = M®P, (102) 

where all the operators are understood to be restricted to the subspaces spanned by |§ ± (x)). 
Obviously, this commutator expression is satisfied whenever the following two conditions are 
separately satisfied: 

[R((p),M]=0, [R((p),P]=0, (103) 

where the arbitrariness of the angle q> has been exploited. But, again, this condition is sufficient 
but not at all necessary. In fact, it is not difficult to see that if one choose either 
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(105) 

with mi,m2,pi i p2 € C, then Eq. (1100b is always satisfied but Eqs. ( 1103t are satisfied only by 
the first choice ( 11041 ) that represent orthogonal matrices whenever m\ + m\ = 1 = p\ + p\. 
It is worth to notice that a circular polarizer (left or right) is represented by the Jones matrix 



1/1 ±i 
2\Ti 1 



(106) 



which is of the form ( 11041 i. Conversely, a half-wave plate with fast axis at angle a with respect 
to the horizontal axis is represented by the Jones matrix 



cos(2a) sin(2a) 
sin(2a) — cos(2a) 



(107) 



which is of the form (1105b . 



